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ABSTRACT 

Ultraviolet spectropolarimetry of three bright high-redshift low polarization 
quasars (LPQs) was obtained with the Faint Object Spectrograph of the 
Hubble Space Telescope. Two of the quasars, PG 1634+706 and PG 2302+029, 
had polarizations p ~ 0.5 — 1.0% throughout the ultraviolet, and showed no 
significant variation of polarization amplitude or position angle with wavelength. 
PG 2302+029 was also marginally (2.4cr) circularly polarized in the optical 
continuum. For the highest redshift quasar, PG 1222+228 (Ton 1530), the 
polarization was measured down to rest wavelengths below 800 A. Although 
the continuum of PG 1222+228 was weakened by Lyman limit absorption from 
an intergalactic gas cloud, the polarization increased sharply from 1% to about 
4.5%, a change of 4(7 significance. This abrupt rise in polarization does not 
appear attributable to any known instrumental artifact. 

These UV polarizations were only slightly less than those previously observed 
for these same objects in the optical. The polarization spectra were flat with 
a typical slope of the polarized flux pF„ oc i/ -0 - 8:t0 - 5 . Unlike the polarization 
spectra of several high luminosity Seyfert 1 nuclei studied previously, such a flat 
wavelength dependence of polarization cannot be explained by scattering from 
dust grains. The hypotheses that the polarization in these quasars is produced 
by transmission through aligned interstellar grains (in the Milky Way or the 
host galaxy), or by a synchrotron power law component, also appear to be 
ruled out. These observed spectra are consistent with a wavelength-independent 
polarization proportional to the total nonstellar light or, possibly, to the 
contribution of the blue thermal component. 

The polarization spectra have insufficient signal-to-noise to locate the 
scatterers with respect to the continuum source and the much larger broad 
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line region. A decrease in amplitude and rotation of the position angle of the 
polarization vector at the shortest wavelengths, which could result from General 
Relativistic effects near a spinning black hole, was not observed. In fact, in PG 
1222+228, the polarization was observed to increase at the shortest wavelengths. 
The rise in polarization with frequency is so sharp that it cannot be due to 
any wavelength-independent polarizing mechanism at any radius in an accretion 
disk. Such a rise could be attributable, for example, to a relative increase in 
scattering opacity over absorption at higher frequencies. 


Subject headings: nuclei — galaxies: polarization — quasars 
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1. INTRODUCTION 

Spectropolarimetry is a powerful diagnostic tool for studying the continuum and line 
emitting regions of the variety of active galactic nuclei (AGNs). The dominant polarization 
observed in blazars from millimeter to ultraviolet wavelengths, for example, is believed 
to originate from a single, beamed, synchrotron component (Landau et al. 1986, Impey 
and Neugebauer 1988, Bregman 1990). In other types of active galaxies, the polarization 
is lower (p < 1%), and its interpretation is more difficult. Among radio galaxies and 
Seyfert galaxies, there is a vigorous debate on the relative importance of polarization from: 
transmission through aligned dust grains, dust scattering, and electron scattering (Miller 
and Goodrich 1990, Cimatti et al. 1993). In many cases the interpretation is complicated 
by the effects of dilution from an underlying host galaxy. The polarization position angle 
is usually independent of wavelength, suggesting that the polarization in a given object 
originates in a single physical process. Depending on the polarization of the emission lines 
and the continuum, the scatterers are believed to be located near one of the following: the 
continuum region (sub-parsec scales), the broad line region (parsec scales), or the narrow 
line region (kiloparsec scales). 

The degree of continuum polarization in normal quasars is generally low, with a weak 
tendency to increase towards shorter wavelengths (Stockman, Moore, and Angel 1984, Mead 
et al. 1990). These low levels of polarization, combined with the faintness of the sources, 
has resulted in little data being published. Webb et al. (1993), hereafter WMSI, made 
broad band polarization measurements through U , B , V, R and / filters of fifteen bright 
quasars. They compared the results with several models for the origin of the polarized flux. 
The different models included synchrotron emission, scattering from electrons or different 
types of dust grains, transmission through aligned grains, and scattering from electrons in 
an accretion flow. In 9 of the 15 quasars studied, the polarization was best understood as 



- 6 - 


scattering from dust grains. These objects tended to be redder and, in several cases (e.g. 
Mkn 486), had broad emission lines which showed the same polarization as the continuum. 
For three objects, the best fit to the data was a model with polarization proportional to 
the flux of the blue thermal component. These latter objects tended to have lower levels of 
polarization and unpolarized (where measured) broad emission lines. 

Differentiating between the different models for the origin of polarization can best be 
accomplished with data spanning a wide range of wavelengths. Polarimetry of high-redshift 
quasars is particularly desirable because the degrees of polarization predicted by the 
different models diverge in the ultraviolet. 
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2. OBSERVATIONS 

Quasars selected for observation with the Hubble Space Telescope (HST) were chosen 
for their high redshifts, brightness, and measurable polarization at optical wavelengths. 
The high redshifts ensured that the continuum polarization could be measured far into the 
rest-frame ultraviolet. All three quasars have strong ultraviolet excesses, and existing high 
accuracy ground-based polarimetry ( a(p ) = 0. 1-0.2%). PG 1222+228 has a radio power in 
excess of 10 26 W Hz -1 ; the other two have been detected at 8 GHz but are in the range of 
radio power generally referred to as radio-quiet (Kellerman et al. 1989). The three quasars 
and their properties are listed in Table 1, along with a summary of the observations. 

Ultraviolet spectropolarimetry of PG 1222+228, PG 1634+706, and PG 2302+029 
was made with the Faint Object Spectrograph (FOS) of the HST on 11-12 Nov 91, 23-24 
Jul 1992, and 28 Jul 1992 (UT), respectively. All three quasars were observed with the 
G270H grating, which gives coverage from about the atmospheric cut-off down to 2200 A. 
PG 1222+228 and 1634+706 were also observed with the G190H grating covering from 
2300 A down to 1600 A. The continuum fluxes of these two quasars are strongly suppressed 
by optically thick intergalactic Lyman limit absorbers at wavelengths below 2320 A and 
1820 A, respectively. PG 2302+029 was not observed with G190H for lack of time. The 
approximate dispersion for G270H is 0.52 A per pixel and for G190H is 0.37 A per pixel. 
All observations used the blue Digicon detector. Total integration times were 14400 s 
(G190H) and 12960 s (G270H) for PG 1222+228, 5280 s (G190H) and 5280 s (G270H) for 
PG 1634+706, and 7920 s (G270H) for PG 2302+029. 

The FOS serves as a polarimeter by rotating a waveplate/Wollaston prism combination 
into the beam just behind the entrance apertures. These observations used the B waveplate 
which is optimized for observations at ultraviolet wavelengths. The polarimeter has been 
described by Allen and Angel (1982), and calibrations are given by Allen and Smith 
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(1992). A polarization observation consists of eight integrations at different positions of the 
waveplate; each polarization observation was completed within a particular orbit. Each 
target was acquired using an FOS binary search in a 4.3 arcsecond square aperture, and 
polarimetry was carried out using the same aperture. While this aperture allows almost all 
of the light from the aberrated (pre-COSTAR) HST images to enter the instrument, only 
55% of the photoelectron image of the source is actually sampled by the diode array in the 
Digicon. Since Allen et al. (1993) have shown that the instrumental polarization of the 
FOS polarimeter is under 0.1%, it will be neglected in what follows. 

The data were reduced using a program written by Richard Allen of the FOS 
Instrument Team. This program was used to correct the data for instrumental effects. A 
correction was made for dead diodes and the spectra were scaled to counts per second per 
pixel. The sky contribution was sufficiently small that sky subtraction was not required; 
however, a low level of instrumental background noise was subtracted. The data were 
flat fielded, then flux calibrated through multiplication by an inverse sensitivity function, 
and then wavelength calibrated. The retardation coefficients, Wollaston angles, and the 
waveplate angles were combined with a least squares fit to the data points at each waveplate 
position to produce the Stokes parameters. 

The data from each of the two Digicon y-steps were then merged and a correction was 
made for the variation with wavelength of the position of the fast axis of the waveplate. 
Lastly, the resultant Stokes parameter values were averaged over broad wavelength bins 
due to the low polarization and the faint signal. Values of linear polarization have been 
corrected in all cases for the bias that results at low signal-to-noise levels (Wardle and 
Kronberg 1974). Tests showed that the resulting error bars were statistically well-behaved, 
in the sense of decreasing with the square root of the number of channels averaged. Field 
star polarimetry is available only for PG 1222+228, from WMSI. Since the interstellar 
polarization close to the line-of-sight is very small ( <0.3%), and its wavelength dependence 
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in the ultraviolet is poorly known, no correction for interstellar polarization is made in this 
paper. 

Although our motivation was spectropolarimetry, these high signal-to-noise observations 
produced high quality ultraviolet spectra of the three quasars. In a subsequent paper, we 
will consider these spectra, which cover the line-rich region of the Lyman-alpha “forest” 
in two quasars. Here we present the HST spectropolarimetry, and combine it with 
ground-based spectroscopy and polarimetry, in an attempt to constrain models for the 
origin of the polarized flux. 


2.1. Spectropolarimetry 

The measured polarizations were low and wavelength-independent, with the exception 
of PG 1222+228 in the extreme ultraviolet. Table 2 gives a summary of the ultraviolet 
polarimetry, with the Stokes parameters and derived values of polarization and position 
angle calculated in broad wavelength bins. Figure la-c shows the total intensity /, and the 
normalized Stokes parameters Q/I and U/I in adjacent vertical panels for each quasar. 
The total intensity is binned by two pixels for each quasar. The light dashed lines show 
Stokes parameters binned by 32 pixels, which corresponds to 16.6 A for the G270H grating 
and 11.5 A for the G190H grating. The filled circles plus error bars show the more heavily 
binned Stokes parameters, which were calculated using the wavelength ranges listed in 
Table 2. 


2.1.1. Line and Continuum Polarization 

An important motivation for the spectropolarimetry of quasars is the ability to 
distinguish between the polarization of the continuum and that of the strong emission lines. 
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In principle, this allows one to locate the polarizing medium with respect to the broad 
line region, as well as with respect to the much smaller continuum source. In practice, 
this separation is difficult due to the low signal-to-noise of the spectropolarimetry on faint 
objects. Results for PG 1634+706 and PG 2302+029 are shown in Table 3, where the line 
polarization is measured at the position of redshifted Lyman-a, and the continuum brackets 
the line. Table 3 also shows a comparison using a more extended region of continuum, 
given in the third line of Table 3 for each object, which gives a continuum polarization of 
better signal-to-noise. For PG 2302+029, this second comparison includes Lyman-a, C III 
and C IV in the line measurement, and the intervening regions were used for the continuum 
measurement. 

The polarization values of the line and the line plus continuum axe similar; and 
are consistent with the hypothesis that line and continuum have the same polarization. 
However, it is equally important to ask whether the hypothesis that the Lyman-a line is 
unpolarized is inconsistent with the data. In PG 1634+706, the Lyman-a equivalent width is 
118 A. If this line were unpolarized, its flux would dilute the polarization in the wavelength 
region of the line down to 0.09%. The line region in fact has a polarization of 0.23 ± 0.21%, 
so the hypothesis that the line is unpolarized cannot be ruled out. In PG 2302+029, the 
Lyman-a equivalent width is 71 A. If that line were unpolarized, its flux would dilute the 
polarization in the wavelength region of the line down to 0.58%. The line region in fact 
has a higher polarization of 0.86 ± 0.31%, but the difference is not significant. If it is 
assumed that all three strong emission lines are unpolarized, the expected polarization of 
the summed line regions is 0.42%. This can be compared with the measured polarization of 
the summed line regions of 0.72 ± 0.23%, which again is not significantly different. No firm 
conclusion on whether or not the lines are polarized can be drawn from this data. A factor 
of two better signal-to-noise would be sufficient to make this determination. 



- 11 - 


2.1.2. Circular Polarization 

Circular polarization has rarely been detected in the optical continuum of any active 
galaxy. A number of quasars show significant circular polarization at radio wavelengths 
(Komesaroff et al. 1984), where it is attributed to helical magnetic fields or anisotropic pitch 
angle distributions in a synchrotron source. Moore and Stockman (1981) set upper limits 
of 0.5% on the optical circular polarization in a number of blazars. Recently, Valtaoja et al. 
(1993) have quoted the marginal detection of circular polarization in two blazars. Over the 
past twenty years, though, very few high precision observations have been made. 

The most sensitive test for circular polarization involves binning the entire spectrum 
to produce one number for a quasar. The results for the G270H observations were: 

V/I = 0.0023 ± 0.0018 (PG 1222+228), V/I = -0.0005 ± 0.0009 (PG 1634+706), and 
V/I = -0.0039 ± 0.0016 (PG 2302+029). Thus there is a marginal 2.4<7 detection for PG 
2302+029, significant at the 98% confidence level if the errors are normally distributed. 

Instrumental circular polarization in the FOS is no larger than 0.1% (R. Allen, private 
communication) and the only way that spurious circular polarization can be generated in 
the data reduction is if the signal level is very low, which is not the case for any of these 
quasars. Splitting the spectrum into the same line and continuum regions described in the 
section above gives the following result. Circular polarization in the continuum is -0.45 ± 
0.18%, but circular polarization in the wavelength region of the three strong lines is only 
-0.05 ± 0.18%. If this low level of circular polarization is real, it may be associated with 
the continuum of PG 2302+029. 


2.1.3. Is the Polarization Variable? 
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It is important to determine whether or not the low level of polarization in these 
quasars is variable. Synchrotron polarization in strong radio sources is usually highly 
variable (Angel et al. 1978) because the emissivity reflects nonequilibrium conditions in a 
compact region of a relativistic jet (Blandford and Konigl 1979). By contrast, if polarization 
is caused by scatterers located beyond the broad line region, variations should have low 
amplitude and long timescales. 

Berriman et al. (1990) published white light polarization measurements of PG 
1222+228. The values were p = 0.84 ± 0.24%, 9 = 150 ± 8° (1978), and p = 1.07 ± 0.28%, 

9 = 143 ± 8° (1979). WMSI published an R band measurement from March 1984 of 
p = 0.82 ± 0.20%, 8 = 147 ± 7°. There is no indication of variability in these three 
measurements. PG 1634+706 has not exceeded 0.5% in the two measurements from 
Berriman et al. (1990) and this paper, and the two measurements of PG 2302+029 from 
1983 do not differ significantly. The data are sparse, but there is no evidence of the 
variability that is a hallmark of a nonthermal process for producing polarization. 

2.2. Notes on Individual Quasars 

2.2.1. PG 1222+228 (Ton 1580) 

The HST flux levels in the total intensity spectrum are systematically lower by about 
30% than those measured with IUE. Since the absolute flux calibration of the latter is more 
reliable, we have used the IUE points in our fitting. The highest frequency points drop 
below the level of the true continuum due to the merging of optically thick Lyman-limit 
absorbers below a rest wavelength of 912 A. There is no evidence for near infrared variability 
in 5 observations made over a period of 7 years Neugebauer et al. (1989). 

The most remarkable feature of the spectropolarimetry is the sharp rise in the degree 
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of polarization of PG 1222+228 (Ton 1530) at the shortest wavelengths. As mentioned 
previously, instrumental polarization in the FOS is negligible at all ultraviolet wavelengths. 
Due to the low continuum level, it was necessary to verify that an appropriate value for 
the background was used in the reduction process, as the level of the background has 
a significant effect on the resultant values of linear polarization. The background must 
be subtracted off before Stokes parameters are fit to the modulated intensity measured 
through different orientations of the waveplate. If the background is either overestimated 
or underestimated, so is the polarization. The polarization in the bin from 1994 A to 2319 
A was therefore calculated with three different assumptions for the background level. Since 
the background is electronic rather than astronomical, these experiments were carried out 
on the raw data in units of counts per pixel, and the results were converted into intensity 
units using the sensitivity function. 

The FOS background is due to dark current in the Digicon detectors. It is measured 
periodically and shows little variation, except during passage through the South Atlantic 
Anomaly (when no science observations are made). The validity of the standard background 
was tested by considering the values of intensity, linear polarization, and polarization angle 
which were calculated after the subtraction of three different levels of background. The first 
was the standard value, which corresponds to a count rate of 0.0078 counts s -1 pixel -1 . As 
a minimum, zero background was assumed. As a maximum, it was assumed that all the 
signal in the bluest quarter of the G190H spectrum was background, corresponding to a 
count rate 0.011 counts s -1 pixel -1 . This is likely to be too high for a background level, 
since there is no evidence that the dark current varies by this much during normal science 
observations with the FOS. It is also likely that some of the signal from 1600-1800 A is real 
continuum from the quasar. As evidence for this, the average intensity in the bluest quarter 
of the G190H spectrum of PG 1634+706, where the continuum is extinguished due to an 
optically thick intervening absorber, is only 0.0025 counts s -1 pixel -1 . 
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The results of this calculation are shown in Table 4. The systematic error on 
the derived polarization is given by the range under the assumption of minimum and 
maximum background levels. Even under the unrealistic assumption of zero background, 
the polarization only reduces to 2.7%. It was also possible to test for the rapid increase in 
polarization with the blue end of the data taken with the G270H grating. The polarization 
in the wavelength range 2381-2644 A, redward of the Lyman limit system, is 0.85 ± 0.47%. 
The polarization in the narrow window 2224-2286 A, blueward of the Lyman limit system, 
is 9.23 ± 2.64%. High ultraviolet polarization is seen in observations with two different 
gratings, and it is not an artifact of the assumed background level. The high extreme 
ultraviolet polarization of Ton 1530 appears to be real. 

PG 1222+228 was by far the highest redshift object in the sample, and similar increases 
in polarization at equivalent rest-wavelengths in the other two objects are not ruled out. 
The rise in polarization occurs near the position of a Lyman limit absorber at z = 1.486. 
There are absorbers at other wavelengths (where there is no increase in polarization), and 
the occurence of the Lyman limit absorber at the wavelength where the polarization turns 
up may be a coincidence. 


2.2.2. PG 1634 + 706 

Since this quasar is nearly an order of magnitude brighter than the others, it has the 
most accurate polarimetry. This advantage is offset by the very low intrinsic polarization - 
generally less than 0.4-0. 5% at all wavelengths. 

The total intensity flux spectrum is in good agreement with spectra from IUE, except 
that the higher SNR shows a richly populated Lyman-a “forest” blueward of the Lyman-a 
emission line. The 100 A region of overlap between the G270L and G190H spectra shows 
excellent agreement. This quasar shows no evidence of variability at 2^m (Neugebauer et 



- 15 - 


al. 1989). Optical spectrophotometry, from Neugebauer et al. (1987) was not obtained 
simultaneously with the UV data, and shows a level ~ 20% higher. Rather than artificially 
adjusting any data, in our model-fitting we have interpolated over this small discontinuity, 
which produces most of the residuals in the fit. 

2.2.3. PG 2302+029 

The HST total intensity spectrum for this quasar shows the worst disagreement with a 
previous IUE measurement, which is roughly a constant 30% higher. Evidence that this is 
due to real variability comes from optical observations by Neugebauer et al. (1987), which 
show a level 40-50% higher than the optical spectrophotometry we obtained with the Lick 
3-meter telescope in August 1987. Near infrared monitoring shows that this quasar haw a 
more than 80% probability of being variable at 10/zm (Neugebauer et al. (1989)). We have 
used the Lick optical spectrophotometry in our model-fitting, because it agrees well with 
the IUE data, which was also used. 
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3. MODEL-FITTING AND RESULTS 
3.1. Continuum Fits and Component Fluxes 

The predictions of the different models for the origin of polarization require knowledge 
of the different continuum components contributing to the total flux. Therefore, we fit 
component flux models to the multiwavelength (infrared to ultraviolet) continuum energy 
distributions for each of the program quasars, following the procedures of WMSI. Given the 
high redshifts and high nonstellar luminosities, we are confident that the contamination by 
starlight below 3000 A is less than 5% in each case. 

The principal optical/UV component is assumed to be optically thick emission from 
an accretion flow. We described this blue thermal component with the thin disk models 
developed by Malkan (1983) and Sun and Malkan (1989). The fitted disk parameters, 
black hole mass and accretion rate, assuming zero inclination, are listed in Table 5. The 
long-wavelength (A > 1/im) continuum, as would be appropriate for thermal emission, is 
fitted with a power law with high-frequency cut-off around A = 3/im (see Malkan 1991 
and discussion in 4.2 of WMSI). The continuum at rest wavelengths of 2500-3700 A) is 
contaminated by contributions from the Broad Emission Line Region (mostly a blend of 
Balmer continuum and Fe II line emission). An estimate of this component is included in 
our fits. Figures 2a-c show the continuum fits with a thermal infrared continuum. Table 6 
gives the flux by component in each waveband. 

3.2. Polarization Model Fitting 

Each of the quasar’s observed polarization spectra were fit with seven possible models 
for the origin of the polarization: a power law synchrotron emitter, three types of scattering 
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due to dust, a non-spherical distribution of electrons scattering all of the nuclear light, 
electron scattering of the blue thermal component, and transmission through grains in the 
host galaxy. These polarized flux fits are illustrated in Figure 3a-c. The results of a % 2 
calculation showing which models best fit the data are given in Table 7. 

Our new results confirm several conclusions made by WMSI about quasar polarization. 
The absence of any wavelength-dependence of polarization position angles is consistent 
with WMSI’s findings. In no case can the polarizations be produced by transmission 
through aligned dust grains (the “IT” hypothesis). The peaked polarization described by 
the Serkowski law does not match the observations of any quasar. Nor are any of our 
quasar polarizations well fitted by polarization from small dust grains, with scattering 
cross sections proportional to A -1 or A -2 . There is also poor agreement with the predicted 
polarization produced by dust scattering obeying the “MRN” relation used by WMSI to 
describe normal interstellar grains (Mathis et al. 1977). 

Another possible source of polarization is a nonthermal power law (e.g., synchrotron 
emission). However, nonthermal emission could not produce much of the observed flux at 
all wavelengths, because our multiwavelength continuum fits do not reproduce a power law 
with such a flat slope. Also, there is no evidence for the large amplitude variability that is 
usually associated with blazar emission (Neugebauer et al. (1989)). The possible exception 
is PG 1634+706. Its spectrum in the UV range we observed happens to be rather steep 
(F„ oc v~ 1A )] therefore, the constancy of observed polarization with wavelength results in 
an unusually steep drop in polarized flux with wavelength, which could be consistent with 
the presence of a weak blazar (as WMSI found for PG 1048—09 = 3C 246). This steep rise 
of polarization to longer wavelengths in PG 1634+706, a radio-weak object, could probably 
be ruled out by optical or infrared polarimetry. 

Excellent fits to the data are provided for all objects by the “Nonstellar”, “Nuclear”, 
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and “Thermal” models. These favored explanations for the quasar polarizations all require 
some kind of electron scattering of the nuclear continuum, which we discuss next. 


4. DISCUSSION 

The WMSI broad-band optical polarimetry study narrowed the permitted models for 
the origin of polarization in AGNs to essentially two: dust scattering or electron scattering 
of the nonstellar nuclear continuum. Since the UV polarization is either flat or declining 
with frequency in all three quasars, the dust scattering models, which predict increased 
polarization with frequency, do not in general give good fits to the data. 

4.1. Interpreting the Polarized Flux Spectrum 

Since the polarizations from optical to UV wavelengths tend to be flat, the slope of 
the polarized flux spectrum is very similar to the slope of the total intensity spectrum - an 
average slope is p u x F u = k x I n view of the wavelength-independent position 

angle, it is simplest to attribute this polarized flux to a single continuum component. One 
possibility is that most of the continuum we observe in the UV comes from a single physical 
component (but not a power law), which has a low, approximately wavelength-independent 
polarization at wavelengths longer than 800 A). Indeed, WMSI fitted the rest-UV (A < 
3000 A) continuum primarily with a single component which has the spectral shape of 
optically thick thermal emission. Since this thermal continuum is the dominant continuum 
component in the UV flux, it may not be surprising that the “Thermal” hypothesis gives 
some of the smallest x 2 values. However, this model only gives a good fit because the 
polarization is independent of wavelength. 



- 19 - 


The success of the “Thermal” hypothesis fits does not prove that the UV continuum 
is thermal emission, nor does it prove it arises in an accretion disk. It does prove that 
if the UV is disk light, it has a very small, nearly wavelength-independent polarization 
(of course, this is true for any single component causing the polarization). Yet another 
possibility, the “Nonstellar” hypothesis, is that all of the nuclear continuum is uniformly 
polarized (perhaps by electron scattering). This hypothesis is only distinguishable from the 
“Thermal” hypothesis at wavelengths where the thermal component is not the dominant 
flux contributor - i.e. , in the red. This is why the model curves predicting polarized flux 
from the “Nonstellar” and “Thermal” models overlap throughout the rest-UV wavelengths 
which we observed. The predictions of these two models deviate in the observed infrared, 
where a standard thin accretion disk would no longer produce all of the observed flux. The 
nonstellar hypothesis predicts greater polarization in the infrared than does the thermal 

model. 

The only difference between the “Nuclear” and “Thermal” models is that the former 
has unpolarized BLR emission. This only has a signficant effect at rest wavelengths around 
3800-2300A, which is too small to make a significant difference to the x 2 values in Table 7. 

4.2. Interpreting the Far-UV Polarization of Ton 1530 

Figure 4 shows the wavelength dependent polarization of 1222+228 (Ton 1530), where 
the data in the overlap region of the G190H and and G270H gratings has been combined 
into a weighted average. The sharp rise in polarization is clear when the data is compared 
to three simple models, with polazation independent of wavelength, increasing as A" 1 and 
increasing as A~ 4 . If the far-UV continuum is emitted from a standard thin accretion 
disk, with surface temperature proportional to r~ 3 /\ higher-frequency photons tend to 
originate from smaller radii. Nonetheless due to the finite width of the Planck function, 
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the total disk flux at a given frequency must have significant contributions from a range of 
annuli. The contributions of these annuli to the total flux at nearby frequencies cannot be 
dramatically different. To take an extreme case, suppose that the increase in polarization 
at shorter wavelengths were due to the increasingly dominant flux contribution made by 
the (highly polarized) inner region of an accretion disk. Could the observed sudden rise in 
polarization be explained by an inner disk region with wavelength-independent polarization? 
A calculation shows that the observed change - polarization at least doubles while frequency 
changes by no more than 20% - is too sudden. Consider a simplified disk model in which 
effective temperature goes as r _3//4 . Then the integrated disk flux out to a radius where 
the temperature drops to Tb is proportional to the integral ^ where the upper 

boundary of the integral, Xb , is r . 

The best possible case would be when the polarization in the disk outside of this 
boundary drops immediately to zero. We must then ask whether or not the value of the 
definite integral can increase by a factor of two (or more) when Xb increases by only 20% 
(corresponding to a 20% increase in observing frequency)? A numerical integration of this 
integral over a wide range of values of Xb shows that the answer is always no. The largest 
percent changes in the integral occur for small values of Xb , where we can approximate 
the integral by (3/5) Even in this most favorable case, a 20% increase in Xb only 
increases the integral by 36%, far short of the 100% increase required. Relativistic effects 
in a more sophisticated accretion disk model do not appear to be large enough to alter this 
result. We therefore conclude that if the observed rise in polarized far-UV flux originates in 
an accretion disk, the region where it is produced cannot generate wavelength-independent 
polarization. This inner region of the disk would have to produce polarization which 
increases with frequency. To produce this intrinsic increase in polarization would probably 
require that some scattering process, such as Thomson or Compton scattering, becomes 
dominant at shorter wavelengths. 
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4.3. Possible General Relativistic Efffects 

By symmetry, the intrinsic polarization from an accretion disk should be either 
perpendicular or parallel to the disk’s polar rotation axis projected onto the sky; the details 
depend on the temperature and surface structure of the disk (Webb and Malkan 1986, 
Coleman and Shields 1990). It has also long been realized (Conners, Piran, and Stark 
1980) that if polarized emission arises near a black hole, the electric field position vector 
will be rotated as the light propagates through the strong gravitational field. This general 
relativistic effect will only affect photons emitted from within a few tens of gravitational 
radii. If higher frequency continuum photons originate closer to the black hole, the 
polarization should rotate with frequency, especially if the hole is spinning rapidly. 

Our spectropolarimetry sets an empirical limit of 6° on the rotation of the polarization 
vector. If we adopt the thin disk models fitted in Section 4, this limit is significant, especially 
for PG 1222+228 which has polarimetry extending all the way down to rest wavelengths 
of 700-800 A. We use the calculations of Laor, Netzer, and Piran (1990) to transform our 
thin accretion disk fits to the average viewing angle of cos i = 0.5 (a disk viewed face-on 
would not be expected to show any polarization), enabling a direct comparison with their 
Figures 1 through 4. As described in equations 2 and 3 of Sun and Malkan (1989), for the 
rapidly spinning (Kerr) hole our fitted black hole masses (which assumed a face-on disk) 
must be increased by a factor of 4, while the physical accretion rate hardly changes. To 
transform to a non-rotating (Schwarzschild) hole viewed at cost = 0.5, our tabulated black 
hole masses must be tripled, and the mass accretion rates must be increased by a factor of 
4. Note that for cos i = 0.5, the inferred accretion rates, for either Kerr or Schwarzschild 
holes, are 11-17% of the Eddington limit for these quasars. 

If we first assume a Schwarzschild hole, the predicted rotations of the polarization 
vector in Ton 1530 are 6° in the shortest wavelength bin measured (an observed wavelength 
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of 2156 A), and 4° in the near-UV bin centered on 3029 A. These small rotations are 
consistent with our data, within the observational uncertainties. However, for a rapidly 
rotating hole, the predicted polarization position angle changes are 22° at 3029 A and 32° at 
2156 A, which are strongly ruled out by the data. The position angles in PG 1634+706 have 
large measurement errors, but in PG 2302+029 the position angle is measured accurately 
down to a rest wavelength of 1290 A. Even to this redder limit, the predicted polarization 
rotation of 15 degrees for the extreme Kerr black hole can be ruled out. The smaller position 
angle rotation of 4° predicted for a Schwarzschild hole cannot be ruled out. Although these 
predictions were made for an accretion disk with a local blackbody spectrum, Laor, Netzer, 
and Piran (1990) predict almost as large changes in position angle when they replace the 
local blackbody spectra with simple atmosphere models. 

The conclusion that the Kerr models predict a larger polarization rotation than we 
observe is not sensitive to the assumed inclination of the disks. Roughly speaking this is 
because the radius from which the observed UV photons arise, in units of the gravitational 
radius, increases only mildly as the disk inclination increases. In edge-on disks, the far-UV 
continuum comes from a larger number of gravitational radii (i.e., where the gravitational 
potential is weaker). However this effect is mostly offset by the fact that the relativistic 
polarization rotations are enhanced for a nearly equatorial observer. The observed limits 
on position angle rotations at short rest wavelengths in Ton 1530 (and to a lesser extent 
PG 2302+029) permit interesting speculations within the context of the standard model of 
accretion onto a supermassive black hole. Either the far-UV continuum arises from more 
than a few tens of gravitational radii, or the black hole is not spinning rapidly, or the black 
hole is less massive than predicted by the standard model. 
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5. CONCLUSIONS 

None of the three observed quasars showed wavelength dependence of the strength 
or position angle of the polarization from optical wavelengths down to a rest wavelength 
of 800 A (over two octaves in frequency). This behavior strongly suggests an electron 
scattering origin for the low level polarization, but it is not clear whether these scatterers 
are intrinsic to the continuum-emitting region or well outside it. The signal-to-noise of the 
measurements is insufficient for us to determine if the broad emission lines are polarized or 
not. There is some evidence for circular polarization of the continuum in PG 2302+029, 
although it is only significant at the 98% confidence level. 

The levels of ultraviolet polarization are similar to or slightly lower than the previously 
measured levels of optical polarization. There is no evidence for variable polarization, such 
as would be expected from a compact synchrotron source. The polarized flux spectra are 
flat with a typical slope pF u oc v ~° ' 8±0,5 . Unlike the polarization spectra of high luminosity 
Seyfert 1 nuclei studied previously, such a flat wavelength dependence of polarization 
cannot be explained by scattering from dust grains. The hypothesis that the polarization in 
these quasars is produced by transmission through aligned interstellar grains in the Milky 
Way or the host galaxy also appears to be ruled out. The spectra are consistent with a 
wavelength-independent polarization proportional to the total nonstellar light or, possibly, 
to the contribution of the blue thermal component. 

An increase in UV polarization shortward of a rest wavelength of 800 A has been 
observed in PG 1222+228 at the 4 a significance level. This sharp rise is not attributable 
to any known instrumental artifact. Although the rise occurs near an opacity edge due 
to an intervening Lyman limit absorber, this is most likely a coincidence. Assuming that 
the far UV emission comes from an accretion disk, the rapid polarization rise cannot be 
produced by any intrinsic process that is wavelength independent. Even if the inner regions 
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of the accretion disk are highly polarized, measurements at any UV wavelength include 
emission from a range of annuli, which limits the slope of the wavelength dependence of 
polarization. A scattering process other than electron scattering must become important at 
high frequencies. 

Assuming the standard paradigm of an accretion disk surrounding a supermassive 
black hole, observations of PG 1222+228 at short rest wavelengths allow a search for an 
effect predicted by strong field General Relativity. Taking reasonable estimates of black 
hole mass and disk inclination, the expected position angle rotation of 20 degrees at the 
highest frequencies is not observed. Such a result could be explained by a black hole which 
is not rotating rapidly or is less massive than predicted by the standard model, or it could 
be explained by high frequency emission which originates from more than a few tens of 
gravitational radii from the black hole. 


We are grateful to Rich Allen, for diligent and patient assistance with the polarization 
reduction software, and to Paul Smith for polarimetric insights. Support for this work 
was provided by NASA through GO Program grant numbers 2524 and 3732 from the 
Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Fig. 1 — 

Total intensity I, and the normalized Stokes parameters Qjl and U/I measured for the 
three quasars. The total intensity is binned by 2 pixels, and the dashed histogram shows the 
normalized Stokes parameters binned by 32 pixels. The filled cicles show average values for 
the wavelength ranges listed in Table 2. 


Fig. 2.- 

Spectral energy distributions for the three quasars. Log frequency versus log flux density 
data points are shown with error bars. The multicomponent continuum fit, which includes 
thermal dust as a power law with an exponential cutoff at log v = 14, is shown as a solid 
line. 


Fig. 3.- 

Polarized continuum fluxes for the three quasars. Ultraviolet data points are from Hubble 
Space Telescope observations, optical data includes published U, B ,V, R and I and white 
light (unfiltered) measurements where available, with error bars. Polarization models axe 
shown with different line types and are named with symbols and abbreviations defined in 
the text. 


Fig. 4. — 

The wavelength dependence of the linear polarization of PG 1222+228, using data from 
both G190H and G270H gratings. Simple models, normalized at the mean wavelength and 
polarization of the G270H data, show no wavelength dependence (dotted line), p <x A -1 
(short dashed line), and p oc A -4 (long dashed line). 



Table 1. Details of the Observations 


Properties 

PG 1222+228 

PG 1634+706 

PG 2302+029 

Redshift 

2.046 

1.334 

1.044 

V Magnitude 

15.5 

14.9 

16.0 

Mb* 

-29.9 

-29.6 

-27.9 

Pradio (W/Hz) b 

1.9 x 10 26 

1.6 x 10 25 

1.8 x 10 24 

Date of Observation 

11-12 Nov 91 

23-24 Jul 92 

28 Jul 92 

HST/FOS Gratings 

G270L/G190H 

G270L/G190H 

G270L 

Integration Times (s) 

14400/12960 

5280/5280 

7920 


•Luminosities calculated for Ho = 50 km 1 Mpc 1 } q 0 = 0.5 
b Radio luminosity from Kellerman et al. (1989) 






Table 2. Ultraviolet Polarization 


Object 

<] ^ 

Q/I±c(Q/I) 

U/I ± a{UjI) 

p ± ( 7(p)* 8 ± a(8) Intensity 1 * 

(%) (deg) (erg/ s/ cm 2 / A) 

PG 1222+228 

1994-2319 

0.0354+0.0093 

-0.0305+0.0093 

4.57+0.93 

160+6 

3.92X10" 16 


2224-2492 

0.0031+0.0060 

-0.0171+0.0060 

1.63+0.60 

140+10 

7.14xl0 -16 


2493-2761 

0.0055+0.0045 

-0.0004+0.0045 

< 0.90 

oo 

9.33X10" 16 


2762-3029 

0.0091+0.0037 

-0.0082+0.0037 

1.17+0.37 

159+9 

1.41xl0 -15 


3030-3295 

0.0168+0.0042 

0.0019+0.0042 

1.64+0.42 

3+7 

1.69xl0" 15 


2224-2761 

0.0045+0.0036 

-0.0076+0.0036 

0.81+0.36 

150+12 

8.23X10" 16 


2762-3295 

0.0133+0.0028 

-0.0027+0.0028 

1.33+0.28 

174+6 

1.55xl0" 16 

PG 1634+706 

1900-2320 

0.0004+0.0021 

0.0029+0.0021 

< 0.42 

~41 

1.17X10” 14 


2224-2492 

0.0034+0.0020 

0.0002+0.0020 

< 0.40 

~2 

1.36X10" 14 


2493-2761 

0.0049+0.0018 

0.0009+0.0018 

0.46+0.18 

5+10 

1.34X10" 14 


2762-3029 

0.0024+0.0016 

-0.0017+0.0016 

< 0.32 

~163 

1.90X10" 14 


3030-3295 

-0.0022+0.0024 

0.0034+0.0024 

< 0.48 

~61 

1.19xl0 -14 


2224-2761 

0.0041+0.0014 

0.0006+0.0014 

0.39+0.14 

4+9 

1.35X10" 14 


2762-3295 

0.0007+0.0013 

0.0003+0.0013 

< 0.26 

~11 

1.55xl0 -14 

PG 2302+029 

2224-2492 

0.0001+0.0033 

-0.0058+0.0033 

< 0.66 

~ 135 

3.53X10" 15 


2493-2761 

-0.0009+0.0029 

-0.0106+0.0029 

1.02+0.29 

133+8 

3.53X10" 15 


2762-3029 

0.0031+0.0032 

-0.0023+0.0032 

< 0.64 

~ 162 

3.09X10" 15 


3030-3295 

-0.0018+0.0039 

-0.0068+0.0039 

< 0.78 

~127 

3.12X10" 15 


2224-2761 

-0.0004+0.0022 

-0.0082+0.0022 

0.79+0.22 

134+8 

3.53x 10 -15 


2762-3295 

0.0006+0.0025 

-0.0046+0.0025 

< 0.50 

~ 139 

3.10X10" 15 


a Polarization values are corrected for low signal-to-noise bias 
b Mean intensity over wavelength range indicated 


Table 3. Line and Continuum Polarization 


Object 

Regions 

Q/I±<j(Q/I) 

U/I ± a{U/I) 

P ± <Kp) 1 
(%) 

9±a(9) 

(deg) 

Intensity 
(erg/ s/ cm 2 / A) 

PG 1634+706 

Line* 

0.0025+0.0021 

-0.0019+0.0021 

0.23+0.21 

161+20 

2.67xl0" 14 


Continuum b 

0.0020+0.0013 

0.0013+0.0013 

0.20+0.13 

16+16 

1.33xl0 -14 


Continuum 0 

0.0019+0.0010 

0.0014+0.0010 

0.21+0.10 

18+12 

1.30x 10 -14 

PG 2302+029 

Line d 

-0.0014+0.0031 

-0.0091+0.0031 

0.86+0.31 

131+10 

4.25x 10 -15 


Continuum 0 

0.0008+0.0029 

-0.0083+0.0029 

0.78+0.29 

138+10 

3.07X10" 15 


Line f 

-0.0008+0.0023 

-0.0075+0.0023 

0.72+0.23 

132+9 

3.72x 10 -15 


Continuum 8 

0.0012+0.0023 

-0.0056+0.0023 

0.52+0.23 

141 + 11 

3.00xl0 -15 

‘Wavelength 
b Wavelength 
c Wavelength 
d Wavelength 
c Wavelength 
f Wavelength 

range 2800-2900 A (Lyman-a) 
ranges 2500-2800 A and 2900-3200 A 
ranges 1900-2800 A and 2900-3200 A 
range 2400-2600 A (Lyman-a) 
ranges 2200-2400 A and 2600-2800 A 

ranges 2400-2600 A, 2800-2900 A and 3050-3200 A (Lyman-a, C III and C IV) 

^Wavelength 

ranges 2200-2400 A, 2600-2800 A 

, 2900-3050 A, and 3200-3300 A 




Table 4. Extreme Ultraviolet Polarization of PG 1222+228 

Background 

Intensity 

P ± <Kp) 

9 ± a{9) 


(erg/s/cm 2 / A) 

(%) 

(deg) 

Maximum 

0.2914+0.0024 xl0~ 15 

5.74+1.25 

158.4+6.1 

Standard 

0.3985+0.0024 xlO' 15 

4.37+0.91 

157.7+5.8 

Zero 

0.6116+0.0024 x 10 -15 

2.73+0.59 

158.4+6.0 




Table 5. Component Fitting Parameters 


Object 

c*PL 

S PL (5400A) 

(mJy) 

S b (3646A) 

(mJy) 

BH Mass 
(M© yr' 1 ) 

Accretion Rate 
(Mq yr- 1 ) 

PG 1222+228 

-1.30 

0.17 

0.02 

6.0 x 10 9 

29.0 

PG 1634+706 

-1.10 

0.50 

0.10 

13.0 x 10 9 

60.0 

PG 2302+029 

-1.60 

0.14 

0.05 

3.2 x 10 9 

11.6 


Table 6. Polarization Hypotheses - Reduced x 2 s 


Object 

Nonstellar 

Nuclear 

Thermal 

MRN 

A" 1 

A -2 

IT 

PG 1222+228 

mm 


1.43 

2.33 

2.15 

18.0 

260 

PG 1634+706 



0.37 

2.11 

4.38 

24.2 

135 

PG 2302+029 

mm 

0.43 

0.37 

1.45 

1.32 

6.49 

50.2 
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